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Abstract 

LOFAR (the Low Frequency Array) is the largest ra- 
dio telescope in the world for observing low frequency 
radio emission from 10 to 240 MHz. In addition to 
its use as an interferometric array, LOFAR is now 
routinely used to detect cosmic ray induced air show- 
ers by their radio emission. The LOFAR core in the 
Netherlands has a higher density of antennas than 
any dedicated cosmic ray experiment in radio. On 
an area of 12 km more than 2300 antennas are in- 
stalled. They measure the radio emission from air 
showers with unprecedented precision and, therefore, 
give the perfect opportunity to disentangle the phys- 
ical processes which cause the radio emission in air 
showers. In parallel to ongoing astronomical obser- 
vations LOFAR is triggered by an array of particle 
detectors to record time-series containing cosmic-ray 
pulses. Cosmic rays have been measured with LO- 
FAR since June 2011. We present the results of the 
first year of data. 

1 Introduction 

The Low Frequency Array (LOFAR) [1 , a distributed 
radio telescope, is taking advantage of the recent de- 
velopments in fast electronics. It is using digital inter- 
ferometry methods to combine signals from a large ar- 
ray of single antennas, rather than a few large dishes, 
to image the northern radio sky. In addition to im- 



aged data, one is also able to use the signals of the an- 
tennas directly in the time domain. After the LOFAR 
Prototype Station (LOPES) successfully detected the 
coherent radio pulses from air showers with energies 
exceeding 10^^ eV [2 , LOFAR is now consequently 
used to continuously measure air showers with an 
unprecedented number of antennas in parallel to on- 
going astronomical observations. 

2 LOFAR 

LOFAR is a digital aperture-synthesis telescope and 
consists currently of 48 stations in five European 
countries. In the core of LOFAR, in the North of the 
Netherlands, there are 24 stations concentrated on 
about 12 km^ with an increasing density of antennas 
towards the middle (see figure fl]). Each station con- 
sists of two fields of high-band antennas (HBA) with 
a frequency range of 110-250 MHz and one field of 96 
low-band antennas (LBA) with a sensitivity between 
10-80 MHz. Both types could be used for cosmic- 
ray observations, but as the HBAs are processed by 
an analogue beamformer before being digitized, using 
them is challenging. 

With the use of ring-buffers of 1 GB for each an- 
tenna (Transient Buffer Boards) the time-series of the 
signals are stored and read-out after a trigger. With 
this feature the detection of cosmic rays becomes fea- 
sible and making LOFAR the densest array for de- 
tection of radio emission from air showers. 




Figure 1: The central part of the LOFAR core. On 
the so cahed Superterp an area of about 0.1 km^ is 
populated with about 600 low-band antennas and 300 
high-band antennas. 



Being a radio telescope, LOFAR undergoes thor- 
ough calibration procedures on astronomical objects, 
from which observations of signals from cosmic rays 
with unknown strength will profit. A shared instru- 
ment facilitates the sharing of measurements, meth- 
ods and expertise. 

2.1 Low-Band Antennas 

The main workhorse for the cosmic-ray observations 
are the low-band antennas, shown in figure [2] The 
band from 10 - 30 MHz is usually filled with radio 
interference, thus there is an optional digital filter for 
this band. Every LBA element is sensitive to two or- 
thogonal polarizations being aligned ±45° to north. 
The two dipoles are held in place by a mechanical 
triangular structure, on top of which a low- noise am- 
plifier is located. Corresponding to the length of the 
dipoles the system is resonant at slightly below 60 
MHz with falling fianks towards the edges of the fre- 
quency band, which automatically suppresses regions 
with higher interference. As the dipoles are not per- 
pendicular to the ground (45°), the antennas are sen- 
sitive to all three polarizations and allow an omnidi- 
rectional monitoring of the sky. 

All LBA stations consist of 96 antennas on an irreg- 
ular grid with increasing density towards the middle. 
For every observation either the inner or the outer 
half or only one polarization of all antennas can be 
used due to signal path limitations. Thus, an antenna 
set has to be chosen when setting up an observation. 








Figure 2: Impression from the LOFAR core. In the 
foreground a particle detector unit is shown. Behind 
it are the dipoles of the low-band antennas from one 
station. The black boxes in the background house 
the high-band antennas. 



2.2 Particle Detector - LORA 

The LOFAR Raboud Air shower Array (LORA) is 
an independent scintillator array used to measure 
cosmic-ray induced air showers at energies from 10^^ 
to 10^^ eV [3 . It consists of 20 detector units, each 
containing two scintillators (0.45 m^, NE 114), which 
are distributed in clusters within the core of LOFAR 
as shown in figure [3J LORA acts as an external trig- 
ger to LOFAR. Whenever a coincidence between a 
number of detectors is measured, a trigger is sent to 
LOFAR. The trigger conditions are adjustable to stay 
below the allowed data rate. 

LORA also provides air-shower parameters in order 
to reconstruct the radio signal [3 . For showers with 
energies of more than 10^^ eV LORA has an angular 
resolution of better than 1 deg^ and an axis position 
uncertainty of about 5 m for events that are within 
the array. The overall energy uncertainty is about 25 
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Figure 3: Schematic layout of the Superterp, the cen- 
tral core of LOFAR. The black diamonds indicate the 
positions of the low-band antennas. The scintillators 
(LORA) are marked by red squares. 



2.3 Data acquisition and observations 

LOFAR is in its main purpose a radio telescope where 
proposal-based astronomical observations are carried 
out. According to the proposal's bandwidth antenna 
set and stations are chosen. As for cosmic-ray detec- 
tion a proposal based observation would deliver too 
little measurement time, cosmic ray observations can 
run in the background in parallel to ongoing obser- 
vations, making use of the Transient Buffer Boards 
(TBBs). In principle, all antennas can be read out via 
TBBs. The limitation is set by the number of TBBs 
and electronics available, which is half the number 
of antennas (48) in a station. There is however no 
limit on the number of stations, as every station is 
equipped with TBBs. 

When the particle detector (or a pulse finding al- 
gorithm) issues a trigger and the ongoing observation 
allows, the TBBs are read out and the data from all 
stations is collected. The search for a radio signal 
from cosmic rays is performed offline. Currently, the 
particle detector is used as default trigger for cosmic- 
ray observations. Whenever four out of five clusters 
of detectors report an event, which is the case if three 
out of four detectors have triggered, a trigger is is- 
sued. Antennas in an area of about 4 km around 
the central core are read out, as it is improbable to 
measure radio signals from showers at larger distances 
that triggered the scintillators. 

Additionally, there are ongoing efforts to develop 
a radio-only trigger. In this approach, pulses are de- 
tected in real-time, by processing the digitized radio 
signal in a FPGA (field-programmable gate array). 




Figure 4: Skyplot of the direction of cosmic rays de- 
tected from July 2011 until June 2012. Full circles 
represent events with a reliable reconstruction from 
the particle detector, events with open squares did 
not pass the quality cut. The direction of the mag- 
netic field is indicated by the blue circle. North is 
90°. 



Currently, a threshold search on the absolute values 
of the voltage ADC samples is implemented, but a 
smarter algorithm is needed to suppress pulses from 
sources other than cosmic rays. Here a large training 
set from data obtained with the scintillator trigger 
will be essential. As reading out all antennas re- 
quires a considerable amount of data bandwidth, a 
radio trigger needs to be highly efficient with respect 
to false triggers. As soon as a radio trigger is oper- 
ational, all LOFAR stations can trigger individually 
during ongoing observations. 

3 Event Processing 

The current pipeline consists of two parts. In a first 
step the individual traces of each polarization are 
cleaned for narrow-band interference and searched for 
short pulses, the latter using the cosmic ray direction 
as reconstructed from the LORA data, if no pluses 
can be found otherwise. If more than ten pulses are 
found, a direction fit is performed (beamforming) to 
better identify also those signals, which are is hid- 
den in the (at LOFAR mostly galactic) noise. Fur- 
thermore, calibration constants are applied to correct 
for time delays within the stations. This step in the 
pipeline is used to identify candidate events on which 
a second step is applied. 

In a second iteration the polarizations are com- 
bined per antenna location, the antenna character- 
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Figure 5: Energy distribution of cosmic rays detected 
with LOFAR from July 2011 - June 2012. The graph 
includes only events, where the energy is reliably 
measured with the scintillator array. 



istics unfolded in combination with a plane-wave fit 
of the direction and the signal projected into polar- 
izations x,y,z. This step will also be the place where 
the absolute gain calibration is applied as soon as its 
results are verified. 
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Figure 6: A signal from a cosmic ray air-shower as 
detected with an antenna at LOFAR in modified 
Analogue-Digital Converter Units (ADU). In addi- 
tion the result of the Hilbert-transformation, the en- 
velope, and the RMS level is shown. 



tings of the scintillator array about one good-quality 
cosmic-ray event is measured per 12 hours of obser- 
vation with the LBAs, i.e. in the order of 200 events 
per year. 



4 Data set 

The current cosmic ray data set is depicted in figure 
m It shows the directions of all events that have been 
identified as candidate events. For further analysis we 
will only use those that also have high-quality shower 
parameters provided by the particle detector, such as 
a core position being contained in the detector array. 
The skyplot shows a clear north south asymmetry 
with a probability of 0.69±0.037 of detecting an event 
from the north in radio emission. As the Earth's 
magnetic field points almost north at LOFAR (89.3° 
azimuth, 67.9° down) this asymmetry is expected if 
the geomagnetic emission mechanism is the dominant 
process [4 . For showers that do not show a pulse the 
probability is 0.479 ± 0.017 to arrive from the north, 
which is not a significant deficit and indicates that 
the trigger is uniform and effective below the energy 
threshold for radio detection. 

The energy distribution of the high-quality events 
is shown in figure |5j Clearly visible is a threshold 
effect towards lower energies and the flux limitations 
of the spectrum at higher energies. The excellent 
conditions at LOFAR with low background ensure 
detections of radio emission from air showers exceed- 
ing 10^^ eV in energy. With the current trigger set- 



5 Air shower properties 

Despite the fact that radio emission of air showers 
has been known to exists since the 1960s O |6] the 
dependencies on shower parameters are still not fully 
understood. However, the theoretical picture is mak- 
ing first steps towards consolidation and is therefore 
in need of high-quality data to support or disprove 
the theoretical concepts [7 . 

With the results from LOFAR, which require a 
minimum of 40 antennas for an event, a better un- 
derstanding of emission mechanisms and the depen- 
dencies is expected. Only when the emission mecha- 
nisms are fully understood, detecting radio emission 
from cosmic rays can be an option for an independent 
large scale detector. 

5.1 Pulse Parameters 

In figure |6] an example of the measured field strength 
as function of time recorded with LOFAR is shown. 
It is obvious that the pulse is unambiguously identi- 
fiable. The agreement of the reconstructed direction 
from the LOFAR data of these pulses with the one 
from the particle detector identifies it as originating 
from the air shower. The trace has been cleaned for 
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Figure 7: Simulated event for the LOFAR core. 
Shown are the expected deviations in time from a 
planar shower front as they would be measured with 
LOFAR. Simulated is a spherical shower front with 
radius of 4 km. The differences can be measured with 
the standard time-resolution of LOFAR. 



radio-interference lines (RFI), which only improves 
the signal-to-noise ratio slightly as the background is 
dominated by the Galaxy at LOFAR. 

These conditions will allow studies of pulse pa- 
rameters such as frequency content and polarization. 
These characteristics are suspected to be very sensi- 
tive to the mass of the cosmic ray and different emis- 
sion mechanisms. Polarization results are expected 
as soon as the unfolding of the antenna pattern has 
been verified. 

5.2 Shower Front 

The shape of the shower front of the radio emission 
is also not fully understood yet. While a plane wave 
seems to be a good first approximation for the re- 
construction of the direction, some experimental ev- 
idence has been hinting at a more complex structure 
[8]. Also from a theoretical point of view a moving 
point-like charge should rather produce a curved or 
conical wave front. Simulations, as for example in 
figure [7) show that LOFAR will be highly suitable to 
identify the structure of the shower front. Already 
sampling with 200 MHz enables a time-resolution 
high enough to see the differences at the edge of 
the shower. Together with a complete interferometric 
phase calibration a resolution of sub-nanoseconds is 
achievable. Again, the high number of antennas will 
enable clear statements on a single-event basis, as the 
high number of antennas close to each other highly 
constrains the shape of the shower front. This will be 



especially relevant, if the shower front is sensitive to 
composition parameters such as X^ax- 

5.3 Lateral Distribution 

The current data set already opens the possibility to 
make statements about the form of the lateral dis- 
tribution of the radio signals from an air shower. A 
typical example event is shown in figure ?? Is shows 
quite nicely that all total signals align in a smooth 
band. The single outlier is due to a miscabled an- 
tenna, which has now been confirmed and fixed in the 
hardware after the observation of this issue in cosmic 
ray data. The spread in signal strength might be re- 
lated to the nature of the emission. Final conclusions 
can however only be drawn after the full calibration 
is in place. A parametrization for the lateral distri- 
bution function (LDF) is underway. 

For now we can conclude that almost all of the de- 
tected events (~ 95 %) show the same type of lateral 
distribution with a flattening within a distance to the 
shower axis of about 100-150 m. Outside of this dis- 
tance almost all events show a steep (exponential) fall 
off. Both features are in principal theoretically pre- 
dicted [9]. Furthermore, the signal strength shows 
a clear dependence on the sine of the angle between 
shower axis and geomagnetic field, as well as a linear 
dependence on the primary energy of the cosmic ray. 

When it comes to more subtle features the conclu- 
sions are only tentative. Many signals in individual 
polarizations (local LOFAR polarizations) show fea- 
tures that hint at the interaction of differently polar- 
ized emission components of the signal. In addition, 
there are first hints at a non-symmetrical behaviour 
of the signal with respect to the shower axis. 

5.4 Core position 

The reconstruction of air shower parameters, such as 
the position of the shower axis, from radio data will 
highly depend on an exact knowledge of the shape of 
the lateral distribution of the signals. But even with- 
out this knowledge the high number of antennas at 
LOFAR enable us to infer a core position. This is es- 
pecially important in cases when the core position is 
due to geometric reasons not or only unreliably avail- 
able from the particle detectors. In figure [??] an ex- 
ample is shown. The core is shifted considerably with 
respect to the core provided by the particle detectors 
and classical methods, such as a weighted mean of the 
radio signals. The only criterion to obtain the new 
core is that the number and size of jumps in the LDF 




240 E 



-100 -50 50 

LOFAR East [meters] 




200 250 300 

Distance to shower axis [m] 



Figure 8: Cosmic-ray event detected by LOFAR. The left plot shows the so called footprint. The circles 
represent the measured total signal (quadratically added amplitudes of all polarizations) in the individual 
antennas, where the size indicates the signal strength. The arrival time of the signal is coded in colour going 
from red (early) to yellow (late). The dark red pentagons represent the signals measured in the particle 
detectors, where size corresponds to signal strength. In the right plot the lateral distribution of the radio 
signal with respect to the shower axis is shown. Kinks in the lateral distribution might be due to non- 
axis-symmetrical behaviour of the signal or due to the not yet fully validated inter-station calibration. The 
energy of the shower as provided by the scintillator array is (2.67 ± 0.80) • 10^^ eV. The inclination of the 
shower is l9 = (40 ±1)°. 
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Figure 9: Geometry for a shower core that is not contained within the array of particle detectors. The red 
cross indicates the geometry as seen by the particle detectors. Here, one clearly cannot estimate the distance 
of the shower core correctly, as the minimization will lead to ambivalent results for larger distances. The 
gray square marks the barycenter (signal weighted mean) of the radio shower, which is largely dominated 
by the highest signals. The blue octagon indicates the geometry from a fit of the radio data only. The right 
side then shows the lateral distribution corresponding to the fitted core positions from the radio data. 



need to be minimized, i.e. the sum of all height dif- 
ferences of two neighbouring signals in distance needs 
to be minimized. This approach is especially useful 
in the early stages of calibration. It gives means to 
test whether jumps in the lateral distribution can be 
removed by confirming an incorrect position of the 
core, rather than a difference in calibration. The ac- 
curacy of this technique can in this early stage not 
be compared with the one from the particle detec- 
tors, but it gives a handle to study also showers that 
are not fully contained within the central core. 

6 Outlook 

With continuously growing statistics the data analy- 
sis will advance. Furthermore, the first stable version 
of full astronomical calibration data is tested and also 
compared to measurements with non-astronomical 
calibration sources. This will lead to a reliable ab- 
solute calibration and the possibility to investigate 
polarization properties of the air showers. 

In addition to a continuation of the data analysis, 
thorough comparisons with simulations of radio emis- 
sions of air showers are planned. First comparisons 
are already underway with CoREAS [10 and EVA 
[TT]. With the high density of antennas LOFAR data 
will be able to put severe constraints on models and 
help to improve the understanding of the emission 
mechanisms. 

7 Conclusions 

LOFAR is the largest array of antennas in the North- 
ern hemisphere. It has been designed to accommo- 
date the observation of cosmic-ray induced air show- 
ers in the radio frequency band parallel to ongoing as- 
tronomical observations. Since June 2011 LOFAR is 
continuously gathering data from air showers, based 
on a trigger by an array of particle detectors. First 
results of air-shower characteristics and signal distri- 
butions indicate that the data taken with LOFAR is 
highly suitable to improve our understanding of emis- 
sion mechanisms and dependencies of the emission on 
air-shower parameters such as the particle type of the 
cosmic ray. 
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